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Abstract The non-Markovian behaviour of open quantum systems interacting 
with a reservoir can often be described in terms of a time-local master 
equation involving a time-dependent generator which is not in Lindblad 
form. A systematic perturbation expansion of the generator is obtained 
either by means of van Kampen's method of ordered cumulants or else 
by use of the Feynman- Vernon influence functional technique. Both ex- 
pansions are demonstrated to yield equivalent expressions for the gen- 
erator in all orders of the system-reservor coupling. Explicit formulae 
are derived for the second and the fourth order generator in terms of 
the influence functional. 
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1. Introduction 

The theory of open quantum systems [1] is concerned with the eUm- 
ination of the degrees of freedom of the environment in order to get an 
equation of motion for the density matrix of the reduced system. For an 
appropriate description of an open system's dynamics that exhibits non- 
Markovian features a primary goal is to derive exact representations for 
the reduced density matrix. A closed expression for the density matrix 
can be obtained for a few analytically solvable models, such as for the 
damped harmonic oscillator and for free Brownian motion [2, 3, 4]. In 
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many interesting cases, however, exact representations for the reduced 
density matrix serve as starting points of a perturbation expansion in 
the system-environment coupUng and of the development of numerical 
integration schemes. 

One possibility of carrying out this program is to derive a time-local 
master equation for the open system's density matrix psit) which takes 
the form 



/C(t) is a time-dependent generator, a super-operator in the reduced 
system's Hilbert space Tis- As can be shown with the help of the 
time-convolutionless (TCL) projection operator technique [5, 6] a mas- 
ter equation of the form (1) indeed exists for small and intermediate 
couplings in the case of factorizing initial conditions. Note that Eq. (1) 
is local in time, i. e. that it does not involve an integration over the 
past history of the reduced system. Due to the explicit time-dependence 
of the TCL generator /C(t), however, it does not lead to a quantum dy- 
namical semigroup and, therefore, the generator need not be in Lindblad 
form. 

A perturbation expansion of the TCL generator /C(i) may be found 
in two different ways. One way is to start from the formal solution of 
the von Neumann equation of the total system and to use van Kampen's 
technique of the ordered cumulant expansion [7, 8]. Another way is to 
invoke the Feynman- Vernon influence functional representation of the 
reduced density matrix [2] and to obtain an expansion of the generator 
directly in terms of the influence phase. Both strategies will be compared 
in this paper and shown explicitly to yield identical expansions of the 
TCL generator. 

To be specific we consider a system-reservoir model in which the reser- 
voir consists of a collection of harmonic oscillators with frequencies Un 
and masses m„. The corresponding coordinates and momenta are de- 
noted by Xn and p„, respectively. The reservoir Hamiltonian is therefore 
given by 



The system-reservoir coupling is described by the interaction picture 
Hamiltonian 



-ps{t)=lC{t)ps{t). 



(1) 





Hi{t) = -aX{t)B{t). 



(3) 



a represents an overall coupling constant, X{t) is an interaction picture 
system operator (not necessarily the position coordinate) and B{t) the 
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interaction picture reservoir variable given by 

B{t) Kn ( Xn COS LOnt + sin COnt ] ■ (4) 

n V mnUJn J 

The constants k„ describe the strength of the coupHng of the reservoir 
mode n to the reduced system. The dynamics of the total system in the 
interaction picture is then determined by the von Neumann equation 

j^pit) = aC{t)p{t), (5) 

where we have introduced the Liouville super-operator C{t). It is defined 
by the relation 

C{t)p = i[X{t)B{t),p\, (6) 

where p is any operator of the combined system. 

The aim is to eliminate the variables of the reservoir to obtain an exact 
representation for the reduced density matrix ps{t) of the open system. 
The starting point is the following formal equation which relates the 
reduced density matrix ps{t) at time t to the density matrix p(Q) of the 
total system at the initial time 0, 

/95(t) =trB|T^exp aJ^dt'C{t') p(0)| , (7) 

where tiB stands for the trace over the degrees of freedom of the envi- 
ronment and T^ denotes the chronological time-ordering operator. 

Let us restrict ourselves here to an initial low-entropy state which is 
given by a product state of the form 



p(0) = ps{0) ® PB. (8) 

Here, ps{0) is the density matrix at the initial time and /j_b is the den- 
sity matrix of the reservoir describing a thermal equilibrium state of 
temperature T which is given by the Gibbs state 

PB = ^exp{-pHB), (9) 

where f3 = l/kBT denotes the temperature and kB is the Boltzmann 
constant. The normalization factor Zb represents the reservoir partition 
function. 
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Time-local master equations and ordered 
cumulants 



We are looking for an appropriate expansion of the TCL generator of 
the master equation (1) with respect to the coupling constant a, 



/C(i) = £a-/C„(t). 



(10) 



n=l 



A general formula for the nth-order contribution JCn{t) to the generator 
of the TCL master equation can be derived by employing a technique 
which was developed by van Kampen for the perturbation expansion 
of stochastic differential equations. To explain briefly this method we 
define for any super-operator TZ of the combined system a corresponding 
super-operator (7^) of the reduced system through the relation 



{n)ps = tTB {nips^PB)}. 



(11) 



The formal representation (7) may thus be written in the equivalent 
form 



°° (T" / 



a 



n=0 



dt'C{t') 
dt'C{t') 



PsiO) 
Ps(0), 



(12) 



where we have expanded the time-ordered exponential in powers of a. 
Explicitly, the first few terms of this expansion take the form 



psit) 



1 + a 



dti(£(ti)) + a' 

■tl /■t2 



dti / dt2{Cih)Cit2)) (13) 



W [ dh r dt2 [ ' dt3{ji:iti)Cit2)Cit3)) 

Jo Jo Jo 



+ ■ 



PsiO). 



Differentiating this equation with respect to time we get 
d_ r rt 

dt' 



-psit) = \a{Cit)) + f dh{Cit)Cih)) 
' L Jo 



(14) 







W / dti / dt2{Cit)Citi)Cit2)) + 







P(0)5. 



The strategy is now to invert the expansion on the right-hand side of 
Eq. (13) with the aim to express ps{Q) in terms of psit), and to substitute 
the result into Eq. (14). As was shown by van Kampen this procedure 
can be carried out in a systematic fashion to yield an expansion for the 
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equation of motion in powers of the coupling. The result is a time-local 
master equation of the form (1), where the nth-order contribution to the 
TCL generator (10) is given by 

ft rtl ftn-2 

JCnit)= dh dt2... dtn-l{jCit)jC{h)C{t2)...C{tn-l)) oc 

Jo Jo Jo 

(15) 

The quantities 

{C{t)C{h)C{t2)...C{tn-l))oc (16) 

^ ^(-i)(^-i) (/:(*) . . . c{u)){cit,) . . . ciu)){c{u) . . . c{tm)){. . .) 

are called ordered cumulants. They are defined by the following rules. 
First, one writes down a string of the form {£...£) with n factors of C 
in between the brackets. Next one partitions the string into an arbitrary 
number of q substrings (1 < g < n) of the form {C . . . C) by inserting 
angular brackets between the Cs, whereby each substring contains at 
least one factor of C The resulting expression is multiplied by a factor 
(-l)(5-i) and all £s are furnished with a time argument in the following 
way. The first factor is always C(t). The remaining Cs carry any permu- 
tation of the time arguments ti,t2, ■ ■ ■ , tn-i with the only restriction that 
the time arguments in each substring must be ordered chronologically. 
In Eq. (16) we thus have 

t>...>ti, tj>...>tk, ti>...>tm, ... (17) 

Finally, the ordered cumulant is obtained by a summation over all pos- 
sible partitions into substrings and over all allowed distributions of the 
time arguments. 

For the thermal state (9) and an interaction which is linear in the 
reservoir coordinates and momenta, as in Eq. (3), we have 

{Cih)Cit2) . . . Cit2n+l)) = 0. (18) 

Thus, only even-order contributions survive. The above rules then lead 
to the following explicit expressions for the second- and the fourth-order 
contributions to the TCL generator: 

lC2{t) = f dti{C{t)C{ti)), (19) 
Jo 

ft ftl ft2 

lCi{t) = dti dt2 / dt3 (20) 
Jo Jo Jo 

{{C{t)C{ti)C{t2)C{h)) - {C{t)C{h)){C{t2)C{h)) 
-{C{t)C{t2)){C{t,)C{U)) - {C{t)C{h)){C{t,)C{t2))) . 
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In general, one expects that a time-local master equation whose gen- 
erator consists of only the first few terms of the expansion provides a 
good description of the reduced dynamics for weak and moderate cou- 
plings. However, it should be emphasized that an expansion of the form 
(10) need not exist for strong couplings. What happens in these cases 
is that Eq. (12) cannot be solved uniquely for psiO). In other words, 
the initial state psiO) is not uniquely determined by the state psi^) at 
time t. Specific examples of the application of this technique to phys- 
ical models and of the breakdown of the TCL expansion in the strong 
coupling regime are discussed in [1]. 

3. Influence functional approach 

In each order of the cumulant expansion the TCL generator involves 
certain combinations of n-point correlation functions of the reservoir 
variables which enter the expressions (16) for the ordered cumulants. 
Another strategy of obtaining an expansion of the TCL generator is to 
eliminate first the reservoir variables completely form the expression (7). 
This is indeed possible for the present model since the initial state (8) 
is Gaussian with respect to the reservoir variables and since the system- 
reservoir interaction is linear in these variables. Following the procedure 
used in [9] one finds the following exact super-operator representation 
of the reduced density matrix. 



Equation (21) provides a complete description of the influence of the 
reservoir on the reduced system. The motion of the system is determined 
by a time-ordered exponential whose complex phase i^t [Xc, Xa] is a 
bilinear functional of the super-operators Xc{t) and Xa{t). The action of 
these super-operators on any matrix p is defined through the commutator 
and the anti-commutator as 



ps{t) = Tt exp (ia'^^t [Xc, Xa]) ps (0) , 



(21) 



where 




xs)p = [x{t),p], 

Xa{t)p = {X{t),p}. 



(23) 
(24) 



T^ denotes the time-ordering of these super-operators Xc{t) and Xa{t). 
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The time-ordered exponential function in Eq. (21) represents the 
super-operator analogue of the Feynman- Vernon influence functional, 
which is usually derived utilizing path-integral techniques [2, 3]. Note 
that the double time-integral in Eq. (22) is already time-ordered for the 
integration is extended over the region t>t'>t"> 0. Two fundamental 
2-point correlation functions enter the above expression for the influence 
phase, namely the commutator function 



D{t-t')=i[B{t),B{t')\, 
and the anti-commutator function 

Di{t-t') = txB{{B{t),B{t')]pB), 



(25) 



(26) 



which are known as dissipation and noise kernel, respectively. 

In order to construct a perturbation expansion for a time-local equa- 
tion of motion for the reduced density matrix wc can proceed in a similar 
way as was done for the cumulant expansion: We first expand the expo- 
nential function in Eq. (21), 



Ps{t) 



E 

m=0 



T^ ($t [X,,Xa]rps{Q) 



a 



1 + ia'^t [Xe, Xa\ - — T^ {^t [Xe, Xa\Y + 



and take the time derivative which leads to the expression: 

d 



dt 



ps{t) 



-m 2m J 



m=l 



(27) 
PsiO), 

(28) 
PS (0) . 



A time-local master equation may again be found by solving Eq. (27) 
for ps{0) within the desired order and by substituting the result into 
the right-hand side of Eq. (28). It is obvious that this procedure yields 
an expansion for a time-local generator of the master equation which is 
identical to the one obtained from the cumulant expansion: The expan- 
sions (12) and (27) are indeed identical to all orders in the coupling and, 
therefore, lead to one and the same expansion of the TCL generator. 
The comparison of the respective expansions shows that the relations 



T^ {^t[Xe,Xa]y 



ml 



z"»(2m)! 



f dt'C{t') 
Jo 



2m\ 



(29) 
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hold for all m = 1, 2, 3, . . . 

Let us illustrate this point by an explicit determination of the sec- 
ond and the fourth order generator. To lowest order Eq. (27) yields 
Ps{t) = Ps(,0). To obtain a second-order equation of motion we thus 
have to substitute this lowest-order expression into the right-hand side 
of Eq. (28) , keeping only the first term. This leads to the second-order 
master equation 

^^psit) = IC2{t)ps{t) (30) 

with the generator 

JC2{t)=ij^<^t[Xc,Xa]= J^dti{Cit)Citi)), (31) 

where we made use of Eq. (29) for m = 1. Obviously, this expression 
coincides with (19). 

The fourth-order contribution of the TCL generator is found by first 
inverting Eq. (27) in second order which gives 

ps{0) = {l-m[Xc,Xa])ps{t), (32) 

and by inserting this expression into the right-hand side of Eq. (28) , keep- 
ing only fourth-order terms. The resulting fourth-order master equation 
is, obviously, 

±ps{t) = [K:2{t)+K:4{t)]ps{t), (33) 
where the fourth-order contribution to the generator reads 

/C4(t) = -^^T^ i^t [Xc, Xa]f + $t [Xc, Xa] j^^t [Xc, Xa] . (34) 

Invoking Eq. (29) for m = 1 and m = 2 this can be transformed into 
)C4{t) = f dh dt2 r dh{C{t)C{ti)C{t2)C{h)) 

Jo Jo Jo 

- f'dh f dt2 dh{C{t)C{ti)){C{t2)C{h)). (35) 
Jo Jo Jo 

Note that the f2-integral in the second term on the right-hand side of this 
equation extends from to t. The triple time-integral of this term may be 
brought into time-ordered form by appropriate substitutions of the time 
variables. More precisely, we may write this integral as a sum of three 
integrals each of which extends over the region t > ti > t2 > t^, > 
In this way, one easily recognizes that Eq. (35) for the fourth order 
generator becomes identical to Eq. (20). 



9 



By use of Eqs. (31) and (34) one easily finds explicit expressions for 
the second and for the fourth order contribution to the TCL generator 
in terms of the dissipation and the noise kernel and of the commutator 
and the anti-commutator super-operators, namely 

]C2{t) = dh |^D(t - h)Xc{t)Xa{h) - ^Di{t - h)X^{t)XSi)^ , 
and 

/C4(i) = ] f dti f dt2 I ' dh 
4 Jo io Jo 

{ [Di{t - t2)Di{ti - ts) + Di{t - t3)Di{ti - t2)] 

xX,{t)X,{h)X,it2)X,{t3) 
- [D{t - t2)D{h - t-s) + D{t - h)D{ti - t2)] 

xX,{t)X,{h)Xa{t2)Xa{h) 

-i [Di{t - t2)Dih - h) + D{t - t3)Di{ti - t2)] 

xX,{t)X,iti)X,it2)Xaih) 

-i [D{t - t2)Di{ti - ia) + Di{t - t3)D{ti - t2)] 

xXc{t)Xc{h)Xa{t2)Xc{t3) 
+D{t - t2)D{h - h)X^{t)Xa{t2)X^{h)Xa{h) 
-Di{t - t2)Di{ti - t3)X^{t)XS2)Xc{ti)X^{U) 
+iD{t - t2)Di{ti - t3)X,(t)Xa{t2)Xc{h)Xc{t3) 
+iDi{t - t2)D{h ~ t-3)X,{t)X,{t2)X,{ti)Xa{t3) 
+D{t - t3)D{h - t2)X,{t)Xa{t3)X,{ti)Xa{t2) 
-Di{t - t3)Di{ti - t2)X,{t)X,{t3)X,{ti)X,{t2) 
+iD{t - t3)Diiti - t2)X^{t)Xa{t3)X^{ti)X^{t2) 
+iDi{t - t3)D{h - t2)Xc{t)Xcit3)Xc{tl)Xa{t2)}. 

4. Conclusion 

The perturbation expansion of the generator of a non-Markovian, 
time-local master equation may be constructed through the technique of 
ordered cumulants or, equivalently, by use of the exact influence func- 
tional expression for the reduced density matrix. We have made explicit 
the relation between both approaches, invoking the connection between 
time-ordered products of the influence phase and of the interaction Li- 
ouville operator. 

It should be clear, however, that an expansion directly through the 
influence functional is only useful, of course, provided an explicit ex- 
pression for the influence functional is available. This was the case in 
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the present study for the reservoir was assumed to be describable by a 
Gaussian (thermal) state. In this respect, the ordered cumulant expan- 
sion is more general since it does not rely on the Gaussian property of 
the environment. 

An important generalization of the present investigation could be to 
include non-factorizing initial conditions for the density matrix of the 
combined system-reservoir state. The corresponding correlations in the 
initial state lead to an inhomogeneity in the master equation. The ex- 
pansion of this inhomogeneity in powers of the system-reservoir coupling 
is known from the time-convolutionless projection operator technique 
[5, 6], while the treatment of non-factorizing initial conditions is also 
possible within the framework of the influence functional technique for 
Gaussian reservoirs [4]. 

Finally we emphasize that the derivation of a time-local generator 
may be important from the numerical point of view. Not only is a 
time-local master equation certainly easier to solve than a generalized 
master equation involving a retarded memory kernel, but it also offers 
the possibility of a stochastic unraveling of the master equation: To all 
orders in the coupling, the form of the TCL generator allows to design 
an appropriate stochastic process {4>{t),ip(t)) for the state vector in a 
doubled Hilbert space TCs © Ti-s such that the average over the quantity 
1^) yields the open system's density matrix [10]. The TCL form of the 
master equation thus gives rise to stochastic wave function algorithms 
for non-Markovian quantum processes. 
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